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ABSTRACT 

Context. The estimated total number of Milky Way globulars is 160±20. The question of whether there are any more undiscovered 
globular clusters in the Milky Way is particularly relevant with advances in near and mid-IR instrumentation. 

Aims. This investigation is a part of a long-term project to search the inner Milky Way for hidden star clusters and to study them in 
detail. GLIMPSE-C02 (G02) is one of these objects, situated near the Galactic plane (7=14° 129, £=-0°644). 

Methods. Our analysis is based on SOFI/NTT JHK S imaging and low resolution (R~1400) spectroscopy of three bright cluster red 
giants in the K atmospheric window. We derived the metal abundance by analysis of these spectra and from the slope of the RGB. 
Results. The cluster is deeply embedded in dust and undergoes a mean reddening of A v ~ 24.8±3 mag. The distance to the object is 
D=4.6±0.7kpc. The metal abundance of G02 is [Fe/H] H 9 6 =-0.33±0.14 and [Fe/H] CG =-0.16±0.12 using different scales. The best fit 
to the radial surface brightness profile with a single-mass Kings model yields a core radius r c =0.70 arcmin (0.9pc), tidal radius r t =15 
arcmin (20 pc), and central concentration c=1.33. 

Conclusions. We demonstrate that G02 is new Milky Way globular cluster, among the most metal rich globular clusters in the Galaxy. 
The object is physically located at the inner edge of the thin disk and the transition region with the bulge, and also falls in the zone of 
the "missing" globulars toward the central region of the Milky Way. 

Key words. Galaxy: globular clusters: general - Galaxy: abundances - stars: distances - stars: abundances 
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1. Introduction 

Star clusters provide us with unique conditions to investigate 
various aspects of stellar astrophysics under tightly "controlled" 
conditions - the clusters are samples of stars with similar ages, 
metallicities and distances. The Galactic globular clusters (GCs) 
can be used to collect information about the formation and early 
evolution of the Milky Way. 

The large area infrared (IR) surveys (i.e. 2MASS, Skrutskie 
et al. l2Q06b have discovered a number of new clusters, hidden by 
the dust extinction in the plane of the Milky Way. These objects 
usually suffer Ay>10-20mag of extinction, making them invisi- 
ble in the optical wavebands. The vast majority of them appear to 
be a few million years old (Ivanov et al. 2002 2005 1 Borissova et 
al. I20031 |2005l 120061 Kurtev et al. 12007] ! but a few have proved 
to be analogues of "classical" globular clusters (Hurt et al. 2000 
Ortola ni, Bica & Barbuy 120001 Kobulnicky et al. 120531 Carraro 
120051 Froebrich, Meusinger & Scholz 2007 ). There are probably 
~10 "missing" globulars in the central region of the Milky Way, 
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based on the asymmetry of the GC distribution (Ivanov, Kurtev 
& Borissova l200"5l 

Spitzer Space Telescope Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire (GLIMPSE, Benjamin et al. |2003l > offers 
an excellent opportunity to carry out an even deeper census of 
such objects than it is possible in the near-IR because of the 
lower extinction at longer wavelengths. A comprehensive search 
for clusters using the point source catalog of GLIMPSE found 
92 candidates (Mercer et al. 2005). As a part of our long-term 
project to find and characterize new Milky Way clusters we stud- 
ied some of them using deep near-IR imaging and low resolution 
IR spectroscopy. In the course of this study we concluded that 
the candidate Nr. 3 in their list (Fig.[TJ is a new Galactic metal 
rich globular cluster. Here we present the evidence and we report 
its properties - metallicity, extinction and distance. We will refer 
to the new cluster as GLIMPSE-C02 (G02), for consistency with 
the designation of Kobulnicky et al. (120051 1 . 

2. Observations and data reduction 

All observations were obtained with Sofl/NTT (Son of ISAAC; 
Moorwood, Cuby & Lidman [T998l on Apr 15, 2006. The instru- 
ment was equipped with a Hawaii HgCdTe 1024x 1024 detector, 
with a pixel scale of 0.288 arcsec px~' . For the spectroscopy we 
used a 1 arcsec slit and the medium-resolution grism, yielding a 
resolution of R~1320 at 2.2 yum. The seeing for all observations 
was 1-1.5 arcsec and the sky was photometric. 
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get spectra), and combined them into final 1 -dimensional spec- 
tra. Finally, we divided them by the spectra we took from the so- 
lar near-analog HIP 59642 (HD 106290) of type Gl V, and multi- 
plied them by a solar spectrum to remove the artificial emission 
lines due to the intrinsic absorption features in the spectra of the 
standard (see Maiolino, Rieke & Ri eke [19961 1 . 




Fig. 1. Pseudo-true color images of G02. It is composed from 
our Sofl/NTT J (blue), H (green) and K s (red) data. The field 
of view is ~3x3 acmin centered at the cluster with coordinates: 
a(2000) = 18 h 18 m 30 s :5 and 5(2000) = -16° 58' 38". North is 
up and East is to the left. 



We collected a total of 16min integration in each of the 
JHKs filters, split into 16 images, jittering within a 3 arcmin 
square box to ensure that there is minimum overlapping of 
the cluster position. Each individual image was the average 
of 3x20 sec frames in /, 6x10 sec frames in H, and 10x6 sec 
frames in K$ ■ The data reduction included flat fielding, sky sub- 
traction, alignment and combination of the individual images. 
The stellar photometry was carried on the final images with 
ALLSTAR in DAOPHOTII (Stetson [19931 . The typical pho- 
tometric errors vary from 0.01 mag for stars with /Ts~10mag 
to O.lOmag for A^lSmag and 0.15mag for Kg ~19 mag. The 
photometric calibration was performed by comparing our instru- 
mental magnitudes with the 2 MASS measurements of about 
1200 stars, covering the color range 0.0<7-/Ts<6.0mag and 
magnitude range 10.0<A's< 15.0 mag. The final photometry list 
contains equatorial coordinates and JHKs magnitudes of 7623 
stars with photometric errors less than 0.15 mag. Artificial star 
tests show that the 80% completeness limit of the photometry is 
at /=18.9 and tfj =17.3 mag. 

The spectra cover the ~2.00-2.35yum region. Fortuitously, 
the slit could be placed to contain 12 stars. The telescope was 
nodded along the slit between the exposures to simultaneously 
observe the targets and clear sky. In total, we obtained 8 im- 
ages of 300 sec. First, we flat fielded them and removed the sky 
emission by subtracting from each field the images from each 
nodding pair. Next, we extracted 1 -dimensional spectra with 
the IRABj task APALL, wavelength calibrated them with the 
NeXe lamp spectra (extracted at the location of the science tar- 



1 IRAF is distributed by the National Optical Astronomy 
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3. Properties of G02 

3.1. Extinction and distance 

The color-magnitude diagram (CMD) of G02 is plotted in Fig. [2] 
The left panel contains all stars with J and Ks photometry. It 
is clearly dominated by field stars. Two main sequences (MS) 
are evident: a nearby, bright one at J—K$ ~0-l mag and Ks ~8.5- 
17.5 mag, and a distant, reddened one at J-Ks~2-4mag and 
Xs~12-17.5mag. There is an indication of a red clump se- 
quence, starting at J—Ks ~3mag, ZTs~13.3mag and extending 
redwards toward fainter magnitudes due to extinction. The clus- 
ter red giant branch (RGB) is distinctly visible at J—K$~4.5- 
5.5 mag in the middle panel where we intentionally selected a 
relatively small region with radius R=60 arcsec near the center 
of G02 to reach good cluster-to-field contrast. Finally, the "field" 
annulus on the right panel, with an inner radius of 138 arcsec 
and an outer radius of 150.48 arsec and identical area, contains a 
mixture of populations and it probably still includes some cluster 
members but we refrain from using a more distant comparison 
area because of the clumpy dust absorption that might severely 
affect the number of field stars and compromise the background 
subtraction. 

The color-color diagram of the cluster field is shown in 
Fig. [3] A comparison of the locus, occupied by cluster stars, with 
the unreddened RGB and MS sequences of Frogel et al. ( 1978) 
yields extinction toward G02 of Ay~24.5+3 mag corresponding 
to a K-band extinction of A^~2.75±0.3 mag. Some of the extinc- 
tion spread might be contributed to differential extinction across 
the face of the clusters, contamination from extended emission 
and background sources. Here and throughout the rest of this 
paper we use the reddening law of Rieke & Lebofsky ( 119851 1. 
Interestingly, the reddening line appears to deviate slightly from 
the sequence of the reddened stars which may indicate anoma- 
lous dust properties in the direction toward G02. 

The cluster luminosity function (LF) is shown in Fig. [3] 
The presence of red giant clump stars is evident and a 
Gausian fit of the bump in the LF gives a mean value of 
/fs = 14.4+0. 15 mag. The clump has an absolute magnitude 
of M K =- 1.83+0.03 mag (Alves [20001 Eqn. 3), yielding a 
distance modulus of (m-M)o~13.45+0.3mag. Here we used 
[Fe/H]=-0.33 (see SectiorTOl). 

The distance to the cluster can also be measured from the 
RGB tip brightness (Ivanov & Borissova 2002). This method is 
hampered by the small number of stars at the upper end of the 
RGB in GCs. Nevertheless, we applied this test as a consistency 
check. Assuming [Fe/H]=-0.33 we obtain an absolute K-band 
magnitude for the tip of M^=-6.8mag (Ivanov & Borissova 
2002). The apparent magnitude of the tip is K$~9.5 mag, yield- 
ing (m-M)o~13.52mag, in excellent agreement with our previ- 
ous estimate. 

We calculated the reddening and distance using the Padova 
isochrones (Girardi et al. 2000) interpolated for metallic - 
ity [Fe/H]=-0.33. The RGB color at the level of the HB 
is J-K-4.15. The true color, obtained from the isochrones, 
is 0.7 which yields E(J-K)-4.05 mag, corresponding to 
Ay=4. 05x6. 2=25.1, consistent with our previous discussion. 
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Fig. 2. Near-IR color-magnitude diagram for all stars on our im- 
ages (left), and the stars in two regions with identical areas - 
a circle with radius R=60 arcsec centered at the cluster center 
(middle), and an annulus centered at the cluster with an inner 
radius of 138 arcsec and an outer radius of 150.48 arsec (right). 
The numbers of stars plotted in each panel are given in brackets. 
The solid line in the middle panel is the best fit to the RGB and 
the crosses mark the stars used to derive the fit in one Monte- 
Carlo realization (see Section |3~2l for details). 



The yielded distance modulus using M^=-1.55 (from Padova 
isochrones) is ~ 12.85, close to both upper values. 

Finally, we obtain the mean values of Ay~24.8±3 mag, 
(ra-M)o~13.3+0.3mag, and D=4.6+0.7kpc for the reddening, 
the true distance modulus, and the distance to the cluster. 

3.2. Metal abundance 

We derived the metal abundance of the cluster using two inde- 
pendent techniques. First, from the spectra of three cluster red 
giants (listed in Table [TJ, we applied the method of Frogel et al. 
(2001), which is based on the behavior of the K-band spectral in- 
dices, IR magnitudes, and colors. The individual reddening and 
true color of each star are obtained by moving the star along 
the reddening vector until it crosses the luminosity class III se- 
quence taken from Frogel et al. dl978l l. All these stars are located 
near the cluster center, and two of them (#1 and #3) lie directly 
on the cluster RGB, making them highly probable members. The 
last one (star #2) suffers from less reddening (AA/j~0.8 mag), 
but the observed differential reddening variations could account 
for this. Indeed, the width of the cluster star locus on the color- 
magnitude and color-color diagram suggests reddening varia- 
tions across the face of the cluster of the order of AAk~0.3 mag. 
Therefore, the reddening of star #2 is identical to that of the other 
two red giants with spectra, within 3<x and differs in less than 
2cr from the mean reddening toward the cluster. Furthermore, 
its angular separation of only 12 arcsec from the cluster center 
strengthens the possibility that this star is a cluster member. 
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Fig. 3. Near-IR color-color diagram for stars in two regions in 
the vicinity of G02, with identical areas - a circle with radius 
R=60 arcsec around the cluster center (open circles; the size is 
proportional to the apparent brightness of the star), and an an- 
nulus around the cluster with inner radius of 138 arcsec and an 
outer radius of 150.48 arsec (crosses). A vector representing vi- 
sual extinction of Ay-25 mag is plotted, using the reddening 
law of Rieke & Lebofsky (119851 1. The face-value colors of RGB 
(solid line) and MS stars (dashed line) from Frogel et al. ( 1978) 
are shown. 

Table 1. Spectroscopic measurements of red giant stars in G02. 



Star 




EW, A 




(J-K)o 


A K 


M K 


[Fe/H] 




Na 


Ca 


CO 


mag 


mag 


mag 




#1 


4.78 


3.62 


15.44 


0.826 


3.002 


-4.828 


-0.29 


#2 


5.88 


2.35 


12.48 


0.601 


2.188 


-4.324 


-0.20 


#3 


4.89 


1.57 


10.33 


0.722 


2.970 


-4.616 


-0.39 



The table contains the equivalent widths of Na, Ca and CO lines, 
true color (J-K)o, individual reddening (A K ), absolute K magnitude 
(M K ), and metallicity in the scale of Zinn. In the calculation of the 
absolute magnitude we used the mean distance modulus to the clus- 
ter: (m-M) ~13.3. 



The spectra are plotted in Fig. [5] We measured the equivalent 
widths of the Na2.20//m, Ca2.26yum, and CO 2.3 fim features, 
and applied the calibration (eqn. 4 of Frogel et al. 200Q that 
takes into account also the stellar color (J-K)o and absolute mag- 
nitude M K . The results are listed in Table Q] Averaging the three 
estimates from the spectroscopy we obtain [Fe/H]=-0.29+0.10, 
in the metallicity scale of Zinn (as implemented in Harris 1996). 
Omitting star #2, we obtain [Fe/H]=-0.34+0. 10. The remaining 
nine stars with spectra were not suitable for metallicity measure- 
ments either because they were too faint or because they were 
early- type field stars. 
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Fig. 4. Normalized near-IR luminosity function, in stars per 
square arc second, for stars within 60 arcsec from the center of 
the cluster (solid line), stars in an annulus with an inner radius 
of 138 arcsec and an outer radius of 150.48 arsec and the same 
area as the cluster region (dotted line), and all stars in the field 
of G02 (long-dashed line). The last function is normalized to the 
whole image area of 7.4 x 7.4 arcmin. 



Next, we derived the abundances from the RGB slope. This 
method was pioneered by Da Costa & Armandroff (1990) and 
it relies on the metallicity-dependent opacities in cool stars. The 
RGB slope is not affected by the reddening, which is an im- 
portant advantage in studies of heavily obscured objects. A crit- 
ical point is to remove the field contamination which can af- 
fect the RGB slope because it represents a population with dif- 
ferent physical parameters. To achieve this we used a Monte- 
Carlo simulation dividing the RGB locus - spanning the ranges 
10.$<K S < 16.8 mag and 4.3<J-K S <5.7 mag - on the CMD into 
bins with sizes 0.5 and 1 .2 mag along the X and Y axis, respec- 
tively. Then, we randomly removed from each bin of the clus- 
ter CMD as many stars as were present in the corresponding 
bin of the field CMD, and determined the cluster RGB slope 
in two iterations, removing the 10cr outliers. This was an inten- 
tionally selected conservative constraint to ensure that we do not 
remove RGB stars or introduce magnitude-related biases. We re- 
peated this process 300 times and calculated the average slope 
and the slope r.m.s.: -0.137+0.012. The calibration of Ivanov & 
Borissova d2002b yields [F e/H]= -0.38±0.14in the scale of Zinn 
(as implemented in Harris |T996l >. and [Fe/H]cG= _ 0.34+0.12 in 
the scale of Carretta & Gratton d!9971 l. 

Given the uncertainties of both methods, the results are con- 
sistent within the errors and we can only state that G02 is among 
the most metal rich GC in the Galaxy. We adopt the average Zinn 
scale metallicity of [Fe/H] =-0.33+0. 14, where the error is tenta- 
tively assigned to the worse error of the two estimates. The cor- 
responding mean metallicity in the scale of Carretta and Gratton 
is [Fe/H] CG =-0.16+0.14. In "The Milk y Way Globular Cluster 
Catalog", 7-th edition of Harris d 19961) one can find only five 
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Fig. 5. K-band spectra of three stars in G02. The position of the 
Na doublet, the Ca triplet, and the first band of CO are indicated. 

clusters with higher metallicity: Lillerl, Terzan5, NGC6528, 
Palomar 10, and NGC 6553. 

3.3. Structural parameters and position in the Milky Way 

To obtain reliable structural parameters and to separate the in- 
trinsic CMD morphology of the cluster from the field we apply 
the statistical decontamination procedure described in Bonatto & 
Bica (2008). Typical cell dimensions used are A/T=0.5mag and 
A(7-//)=A(7-/T)=0.25mag, optimal to allow sufficient star- 
count statistics in the individual cells, without affecting the mor- 
phology of the cluster CMD. The comparison field was the same 
as in the previous sections. 

The structural parameters of G02 were determined using the 
iterative star count method of King ( 1962). The best fit to the ra- 
dial surface brightness profile with a single-mass King's model 
yields a core radius r c =0.70 arcmin (0.9 pc), tidal radius r t = 15 ar- 
cmin (20 pc), and central concentration c=1.33. The radial den- 
sity profile (RDP) of the cluster together with the best fit King- 
like profile are presented in Fig. [6] The bins at 1 .4 and 1 .6 arcmin 
differ from the profile (in one and three standard deviations re- 
spectively) and present an excess of the stellar RPD. The most 
probable explanation of this fact is the variable extinction to- 
ward the cluster caused by the clumpy structure of the absorbing 
matter. Nevertheless, the Kolmogorov-Smirnov tests give a very 
satisfactory probability of 0.97 for this profile fit. 

The cluster is situated very near to the Galactic plane with 
galactic coordinates X=4.46, Y=1.12 and Z=-0.05 in kpc (de- 
fined as in Harris (1996), assuming D=4.6 kpc). It is located 
inside the thin disk, at its inner edge and the transition region 
with the bulge. 

4. Summary 

We found that the object No. 3 from the list of Mercer et al. 
(2005) is a new Milky Way globular cluster, and designate it 
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Radius 

Fig. 6. Stellar RPD for the cluster G02 (open circles). The solid 
curve represents the best fit King-like profile. The radius is in 
arcmin. The X and Y axes are in logarithmic scale. 

GLIMPSE-C02 as the second GC discovered from GLIMPSE 
data. The analysis based on deep JHKs images and moderate 
resolution K-band spectra of three probable members reveals 
a compact, metal rich cluster, with [Fe/H]=-0.33+0.14 in the 
scale of Zinn. The cluster is situated near the Galactic plane be- 
hind Ay~24.8+3mag of visual extinction, inside the thin disk, 
at its inner edge and the transition region with the bulge. The 
mean distance estimated from the red giant clump stars, tip of 
the RGB, and isochrone fits is D~4.6+0.7kpc, placing G02 in 
the zone of the "missing" globulars. 
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